arXiv:1508.04584vl [quant-ph] 19 Aug 2015 


Local optimality of a coherent feedback scheme for 
distributed entanglement generation: the idealized 

infinite bandwidth limit 


Zhan Shi and Hendra I. Nurdin * 
August 20, 2015 


Abstract 

The purpose of this paper is to prove a local optimality property of a recently 
proposed coherent feedback configuration for distributed generation of EPR entangle¬ 
ment using two nondegenerate optical parametric amplifiers (NOPAs) in the idealized 
infinite bandwidth limit. This local optimality is with respect to a class of similar 
coherent feedback configurations but employing different unitary scattering matrices, 
representing different scattering of propagating signals within the network. The infi¬ 
nite bandwidth limit is considered as it significantly simplifies the analysis, allowing 
local optimality criteria to be explicitly verified. Nonetheless, this limit is relevant 
for the finite bandwidth scenario as it provides an accurate approximation to the 
EPR entanglement in the low frequency region where EPR entanglement exists. 


1 Introduction 

Entanglement is a quantum phenomenon in which states (represented by density operators) 
of a composite system composed of several quantum subsystems cannot be written as a con¬ 
vex combination of tensor products of the states of the subsystems. Such entangled states 
have, in recent decades, been of much interest as a resource for quantum information appli¬ 
cations, such as for quantum communication cma- In particular, Einstein-Podolski-Rosen 
(EPR)-like entanglement, generated in the continuous variables such as the amplitude 
and phase quadratures of a Gaussian optical held, has evoked considerable interest over 
discrete-variable entanglement, such as entanglement in finite-level systems like qubits, 
because EPR entangled pairs can be prepared easily and rapidly in quantum optics. In 
this paper, we are interested in EPR entanglement between two propagating continuous¬ 
mode Gaussian fields. Such a kind of entanglement is more accessible compared to EPR 
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entanglement between a pair of single-mode fields produced in, say, inside an optical cavity 

mg. 

EPR entanglement between continuous-mode Gaussian fields can be realized by two¬ 
mode squeezed states produced as the output of a nondegenerate optical parametric ampli¬ 
fier (NOPA). By pumping a strong coherent beam (which can be regarded as an undepleted 
classical light) to a crystal inside the cavity of the NOPA, two vacuum modes of the cav¬ 
ity interact with the pump beam, and photons escaping the cavity through its partially 
transmissive mirrors generate two output beams that are squeezed in amplitude and phase 
quadratures. If the two outgoing fields are squeezed below the quantum shot-noise limit, 
they are considered as EPR entangled beams 0E3- The input /output block representation 
of a NOPA (Gi) is shown as Fig. [lj The NOPA has four ingoing holds and four outgoing 
fields. Among the inputs, £,i OS s,a,i and £,i OS s,b,i are amplification losses, caused by unwanted 
vacuum modes coupled into the cavity. As the two outputs corresponding to the loss fields 
£ioss a i and £,i oss b j are not of interest in this work, they are not shown in the figure. Note 
Fig.0 only presents the ingoing and outgoing noises of interest, and does not show the 
pump beam. 



Figure 1: Input/output block representation of a NOPA. 


In a previous work [6], we have proposed a novel dual-NOPA coherent feedback sys¬ 
tem to produce EPR entangled propagating Gaussian fields, as shown in Fig. [2] It was 
shown that this scheme can produce better EPR entanglement between the propagating 
Gaussian fields £, 0 ut,a, 2 and £, 0 ut,b, 1 (hi the sense of producing more two-mode squeezing be¬ 
tween quadratures of the holds) for the same amount of total pump power used in the two 
NOPAs, and displays more tolerance to transmission losses in the system, as compared to 
a conventional single NOPA and a cascaded two-NOPA system. 

In a subsequent work [7j, we presented a linear quantum system consisting of two 
NOPAs connected to a static passive linear network, realizable by a network of beam split¬ 
ters, mirrors and phase shifters, that are connected in a more general coherent feedback 
configuration, see Fig. [3} Here, the system is ideally lossless, that is, there are no transmis¬ 
sion and amplification losses influencing the system. Hence, each NOPA is simplified to 
have only two ingoing helds, without amplification losses, as shown in Fig. [3} The transfor¬ 
mation implemented by the passive network in this configuration is represented by a 6 x 6 
complex unitary matrix S. By employing a modified steepest descent algorithm, with the 
matrix corresponding to the dual-NOPA coherent feedback network shown in Fig. [2] as a 
starting point, we optimized the EPR entanglement at frequency u> = 0, with respect of 
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Figure 2: The dual-NOPA coherent feedback network. 


the transformation matrix S of the passive network. 
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Figure 3: A coherent-feedback system consisting of two NOPAs and a static passive network 
with six imputs and six outputs from [7]. 


In this paper, we employ the steepest descent method to optimize a coherent feedback 
system shown in Fig. [4j The system contains two NOPAs and a static passive linear network 
described by a 2 x 2 complex unitary matrix S. This system is a more restricted class of 
configuration than the one as shown in Fig. [3j it can be seen that the configuration in 
Fig.g is a special case of the configuration in Fig. [3] Moreover, different from our previous 
work in [7], in which the system shown in Fig. [3 is considered lossless, here we take the 
effect of transmission losses along channels and amplification losses of NOPAs into account. 
However, we neglect time delays in transmission. The effect of delays on EPR entanglement 
generated from related systems can be found in our previous works m- In addition, unlike 
the work in [7j, the system is considered ideally static, that is, we consider the limit where 
the NOPAs are approximated as static devices with an infinite bandwidth [9j. The merits 
of studying this infinite bandwidth limit are twofold: (i) it allows a simplified analysis 
of the system, and (ii) calculations in the infinite bandwidth setting gives a very good 
approximation to the EPR entanglement in the low frequency region, discussed further in 
Section 2.3 In this infinite bandwidth setting, we show explicitly that the choice of the 


scattering matrix in the scheme of |6j is in a certain sense locally optimal with respect to all 
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possible choices of scattering matrices S in the coherent feedback configuration of Fig. [ 4 J 
under certain values of the effective amplitude of the pump laser driving the NOPA. Note 
that there may exist another scattering matrix as a local minimizer that yields better EPR 
entanglement than the network shown in Fig. [2j Searching for such a scattering matrix 
can be a topic for future research. 



Figure 4: A coherent-feedback system consisting of two NOPAs and a static passive network 
with two inputs and two outputs. 


The structure of the rest of this paper is as follows. We begin in Section [2] by giving a 
brief review of linear quantum systems, EPR entanglement between two continuous-mode 
fields, and linear transformations implemented by a NOPA in the infinite bandwidth limit. 
Section [3] describes the system of interest. In Section |4j we discuss the optimization of the 
system. Finally, we draw a short conclusion in Section [5} 


2 Preliminaries 

The notations used in this paper are as follows: 1 = \/ — 1 and Re denotes the real part of 
a complex quantity. The conjugate of a matrix is denoted by • T denotes the transpose 
of a matrix of numbers or operators and •* denotes (i) the complex conjugate of a number, 
(ii) the conjugate transpose of a matrix, as well as (iii) the adjoint of an operator. O mxn 
is an m by n zero matrix (if m = n then we simply write O m ), and I n is an n by n identity 
matrix. Trace operator is denoted by Tr[-] and tensor product is 5(t) denotes the Dirac 
delta function. 

2.1 Linear quantum systems 

Here we consider an open linear quantum system without a scattering process. The linear 
system contains n-bosonic modes a,j(t) (j — 1 ,... ,n) satisfying the commutation relations 
[ai(t),a,j(t)*] = Sij, m-incoming boson fields (i = 1 ,...,m) in the vacuum state, 

which obey the commutation relations [£,i n ,j(t), £mj(s)*] — S(t — s), as well as two outgoing 
fields iout._k{t) (k = 1,2) which are Gaussian continuous-mode fields. A continuous-mode 
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field means that the held contains a continuum of modes in a continuous range of fre¬ 
quencies. Note that, a system may have more than two outputs. However, as we are only 
interested in the entanglement generated by a certain pair of outgoing fields, in this work we 
will only be interested in a particular pair of output holds, labelled out, 1 and out, 2 in the 
following. The time-varying interaction Hamiltonian between the system and its environ¬ 
ment is H int (t) = i(£(t)*L - L*£(t)), where £(f) = [&„,i(£), ■ ■ • L = [L u ..., U} T 

and Lj(j = 1,2, ••• , l) is the j-tli system coupling operator. In the Heisenberg picture, 
time evolutions of a mode aj and an outgoing held operator £ out: i are [UU ITT] : 

=u(ty aj u(t), 

Ut,, (t) =u(ty^(t)u(t), (i) 

where U(t) = ex^ (—i f* H int (s)ds) is a unitary process obeying the quantum white noise 
Schrodinger equation U(t) = -iH- int {t)U(t ). However, this is not an ordinary Schrodinger 
equation as the interaction Hamiltonian H int (t ) is a time-varying observable involving the 
singular quantum white noise processes £(£). This quantum white noise equation has to 
be interpreted correctly within the framework of quantum stochastic calculus, for details 
see mmmm- Employing quantum stochastic calculus, dynamics of a linear quantum 
system is described by quantum Langevin equations and can be written in the following 
form 


z(t) = Az(t) + B£(t), 

£ out (' t ) = Cz(t) + D£(t). 

where 

2 = (af, <) T , 

f — (0 fP e<i ?p ) T 

^iout (s OUt,ll ^>OUt,l^ ^>OUt ,2 7 ^sout^l) 7 

with quadratures mm 

a q j = dj + a*, a P j = (a,j — a*)/i , 

g = b+Cj> $ = (5) 

The linear model described above is ubiquitous in helds such as quantum optics, optorne- 
chanics, and superconducting circuits, and are employed to describe the equations of motion 
for devices as diverse as optical cavities, optical parametric amplihers, optical cavities with 
moving mirrors, cold atomic ensembles, and transmission line resonators, under appropriate 
assumptions on the system’s parameters. 

2.2 EPR entanglement between two continuous-mode fields 

We keep in mind here that in this paper, we investigate EPR entanglement between two 
continuous-mode Gaussian helds rather than entanglement between two single-mode Gaus¬ 
sian helds. In the latter case, the degree of entanglement can be assessed via the logarithmic 


(4) 


( 2 ) 

(3) 
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negativity as an entanglement measure, see, e.g., ra. However, this measure is not di¬ 
rectly applicable to continuous-mode fields. Instead, the EPR entanglement of two freely 
propagating fields containing a continuum of modes, say £ 0 «i,i and £, 0 ut, 2 , can be evaluated 
in the frequency domain by the two-mode squeezing spectra V + (iuj ) and Vh(wu) SHE], 
that will be defined below. 


The Fourier transform of f(t ) is defined as F ( iuj ) = f (t) e lut dt. Similarly, we 

have the Fourier transforms of £, ou t,i(t), £ 0 ^, 2 (t), z(t) and £(f) in ([ 2 ]) and |3| as 5 out) i (tuj), 
Eout ,2 (wu), Z{uS) and S(wu), respectively. Applying (J 2 ]) , (J3J) , we have 


s L,iM + 


sL,iM - 2M 



[10 10] (CZ (icu) + US (?cu)), 



[0 10 — 1 ] (CZ (icu) + D3 (itu )). 


( 6 ) 


The two-mode squeezing spectra V+(uu) and V-(ioj ) are real functions defined via the 
identities 


<(SL,iM + + SL l2 K))) = - a/), 

(7) 

where (•) denotes quantum expectation. As described in [SJ [16], V+(wS) and V-(ioj) are 
easily calculated by, 

V + (iu) = Tr [H 1 (iu)*Hi(xu )], (8) 

v_M= ( 9 ) 

where Hi = [1 0 1 0]H, H 2 — [0 1 0 — 1]H and H is the transfer function 

H(iu) = C(iojI- A)~ 1 B + D. (10) 

The fields ^ out) \ and ^ 0 ut .2 to be EPR-entangled at the frequency co rad/s is [5], 

V(ioj) = V + (iu) + < 4, (11) 

which indicates that the two-mode squeezing level is below the quantum shot-noise limit. 

A perfect Einstein-Podolski-Rosen state is represented by an infinite bandwidth two¬ 
mode squeezing, that is V(icu) = V±(ioj) = 0 for all to. Of course, such an ideal EPR 
correlation cannot be achieved in reality as it would require an infinite amount of energy 
to produce. Thus, we aim to optimize EPR entanglement by making V(iuj) as small as 
possible over a wide frequency range [5j. 

Note that is a sufficient condition for EPR entanglement, with the two beams 
squeezed in amplitude and phase quadratures. However, in general, they may be squeezed 
in other quadratures. Hence, we give the following definition of EPR entanglement. Let 
1 = = e l ^i out)2 with ^ 1,^2 6 (—tt, tt] and denote the corresponding 

two-mode squeezing spectra between and ^tut 2 as V± ^ 2 («^, VhjV^)- 
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Definition 2.1 Fields f ou t,i and font ,2 are EPR entangled at the frequency u> rad/s if 
3 ifl, "02 ^ ( —7T, 7r] sac/i that 

< 4. (12) 

Unless otherwise specified, throughout the paper, EPR entanglement refers to the case with 
V ; i — '02 — 0- EPR entanglement is said to vanish at u> if there are no values of if i and 0 2 
satisfying the above criterion. 


2.3 The nondegenerate optical parametric amplifier (NOPA) 


A NOPA (Gf) is an open linear quantum system containing a two-ended cavity with a pair 
of orthogonally polarized bosonic modes a* and 6 j which satisfy [a;,a*] = S t j, [b^b *] = <5^-, 
[aj, b*\ = 0 and [a*, b 3 \ = 0. By assuming a strong undepleted coherent pump beam onto the 
Xnonlinear crystal inside the cavity, the pump can be treated as a classical field (hence, 
quantum vacuum fluctuations are ignored) and the interaction of the modes a; and b t with 
the pump is modelled by the two-mode squeezing Hamiltonian H — |e (a*b* — aA), where 
e is a real coefficient relating to the effective amplitude of the pump beam, for details see 

unnmnB]. 


As shown in Fig. [TJ interactions between the NOPA and its environment are denoted by 
coupling operators as follows. Modes a* and b t are coupled to ingoing fields fi n , a ,i and fi H} b,i 
via coupling operators L 1 = and L 2 = yfybi, respectively. Unwanted amplification 

losses floss,a,i and fi OS s,b,i impact the NOPA through operators L 3 = y'Tea* and L 4 = 
respectively. The constants 7 and k are damping rates of the outcoupling mirrors (from 
which the output fields emerge from the NOPA), and of the loss channels, respectively. 
Applying Section 2.1 we have the dynamics of the NOPA as jH HH T9], [20] 


di ft) 
h it) 


c-¥) 


0>i ft) T ^ ft) \[Ufin,a,i if) \Z~^floss,a,i ft) , 

hi ft) + “Oj ft) A Jufin,b,i ft) 00 floss,b,i ft) , 


following the boundary conditions HI EH, we have outputs 


£,out,a,i (t) (^) “t“ £,in,a,i (t) 5 

£ out,b,i (t) (t) ^>in,b,i (^) • 


(13) 


(14) 


Define the following quadrature vectors of the NOPA, 


* = [alalblXY, 

p = \p q . p p . p q p p p q pP p q p p ] T 
S Sin,a,2? S in,b,i ’ S loss,a,i ’ S/oss,a,i’ S loss,b,il ’ 

c — \cq cv pq cv 1T 


(15) 
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From ([2j), (j3j) and (10), 


the transfer function of the NOPA is 


H n = 


h 0 

0 hi 

h 2 0 
0 — h 2 


h 2 0 
0 -h 2 

hi 0 
0 hi 


h;i 0 
0 h 3 
h 4 0 
0 —h 4 


h 4 0 

0 —hi 

h 3 0 
0 h 3 


where hj (j = 1, 2, 3,4) are functions of the frequency cu, 


hiiioj) 

h 2 (iu) 

h 3 (iu) 

liiiioj) 


e 2 + y 2 — (k + 2 iu) 2 
e 2 — (7 + k + 2 iu)) 2 ’ 
2 eq 

e 2 — (7 + k + 2 iu>) 2 ’ 
2 Ak / 77 t (7 + k + 2ilu) 
e 2 — (7 + k + 2icu) 2 ’ 
2ey / 7^ 

e 2 — (7 + k + 2 *cu ) 2 


( 16 ) 


(17) 


As reported in HD EH, parameters of the NOPA are set as follows. We set the reference 
value for the transmissivity rate of the mirrors 7 r = 7.2 x 10' Hz. The pump amplitude e 
is adjustable as e = x^ ri where the variable x satisfies 0 < x < 1. We fix the damping rate 


7 = 7 r and set k = Ke with K = 


3xl0 6 


based on the assumption that the value of k 


\/2x0.6X7 r 

is proportional to the absolute value of e and k = 3x ^ when e = 0 . 67 ,.. I 11 this paper, we 
consider the NOPAs have infinite bandwidth case where we take the limit 7 r —)• 00 while 
keeping e and 7 at a fixed ratio ^ — x. In such a case, the transfer function of the NOPA 
in (16) becomes a constant matrix with elements 


(1 - K 2 )x 2 + 1 
x 2 — (1 + Kx) 2 ’ 


x 2 -{l + Kx) 2 ' 
2y/ Kx(l + Kx) 
x 2 — (1 + Kx) 2 ’ 
2x\[Kx 

x 2 - [l+ Kx) 2 ' 


(18) 


It can be seen that for cu e, 7 , k, the constant scalar values of h\ to /r 4 given by (18) 
in the infinite bandwidth limit approximates the frequency dependent values given in (17) 
when the bandwidth is finite. Such an approximation is quite accurate for cu sufficiently 
small, away from e, 7 , k (with no error at cu = 0). Since in practice the EPR entangle¬ 
ment will be in the low frequency region, entanglement in the idealised infinite bandwidth 
scenario provides a good approximation for the entanglement that can be expected in the 
finite bandwidth case. 




















3 The system model 


Consider again the coherent feedback system shown in Fig. [4} The whole network consists 
of two NOPAs and a static passive linear subsystem. The subsystem has two inputs ^' out a 1 
and ti' out ^2 connected to the outgoing fields £, 0 ut,a, 1 of NOPA G\ and £, 0 ut,b, 2 of NOPA G 2, 
respectively. The two outputs £' n b x and £' n a 2 of the subsystem are connected to incoming 
signals °f NOPA G\ and £,i n , a ,2 of NOPA G2, respectively. The incoming fields of the 
system £; niQi i and £^,2 are in the vacuum state m and the EPR entanglement of interest 
is generated between outgoing fields and £out,a, 2 - The transfer function of the passive 

static subsystem is a 2 x 2 complex unitary matrix denoted by S, which satisfies HU 


Sin,a,2 


s 


Son£,a,l 
^>out,b, 2 


(19) 


and 

S*S = SS* = I 2 . (20) 

Also, we shall denote the static passive matrix S corresponding to the dual-NOPA coherent 
feedback network shown in Fig. [ 2 ] as (6J 


S c fb 


0 1 
1 0 


( 21 ) 



Figure 5: Beamsplitter. 


The NOPAs are placed at two distant communicating ends (Alice and Bob). The 
distance between the two ends is d kilometres. Both NOPAs (Gi and G 2 ) have identical 
static transfer functions given by (16) and (18). Transmission loss in each path of the 
network is modelled by a beamsplitter with an unwanted incoming vacuum noise £bSi as 
shown in Fig. [5] The other input is connected to an outgoing held of the NOPAs or the 
subsystem. The outgoing signal £,BS,out of the beamsplitter is the combination of the two 
incoming signals, satisfying C,BS,out = + /3£,bs, where a is the transmission rate and 

/3 is the reflection rate of the beamsplitter, a and /3 are positive real parameters obeying 
0 < a, /3 < 1 and a 2 + /3 2 = 1 [22]. Based on the fact that transmission loss in optical fibre 
is about 0.2 dB per kilometre at telecom wavelengths as reported in [23], the transmission 
rate of each beamsplitter in our system is a = io~°- OOM . 

Technically, transmission losses are accompanied by time delays in the transmission. 
However, here we neglect the time delays in transmission. Nonetheless, the formalism and 
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Heisenberg picture analysis employed here can easily treat the presence of time delays in 
the linear quantum networks considered herein, see [6j [8, [9j, 16] . In previous works on 
related studies [SI [B, US], the effect of time delays has only been to narrow the bandwidth 
over which the EPR entanglement exists, without affecting the EPR entanglement that 
can be achieved in the low frequency region. 

Define the following vectors of quadratures 


z = 

sin 
Sloss 
£,BS = 

f = 

sout 


\CQ cP cl CP ]T 
LS»m,a,l’ Sin,a,15 Sin, b, 2’ Sin,6,2J ? 

\,q £p rq ,p ,p ,q ,p |T 

LSZoss,a,l? S/oss,a,1’ S/oss,6,l? SZoss,6,l? Moss,a,2’ Moss,a,2’ Moss, b, 2’ Moss,6,2J > 

[^BS,o,l> ^BS,a,l ’ ^BS,b,l^BS,b,li ^BS,a,2) &S,2, ^BS,b,2^BS,b,2\ > 

R'? yP cQ cP ] T 
LSowt,l ? Sowt,l 5 Sout,2 1 Sout,2J ’ 


( 22 ) 


Define the real unitary matrix S' as the quadrature form of matrix S'. Based on the 
definitions of the quadratures ([5]), we obtain 


s = \ksk- + 1 k*s*k T , 


(23) 


where 


K = I, 


Note the quadrature form S is, by construction, a unitary sym 
is unitary and symplectic, the latter meaning that S T ^ ^ 


Using the static transfer function of a NOPA given by (16 


(24) 

Hectic matrix. That is, S 
0 1 
-1 0 


S = 


and (18), and given the 


unitary matrix S representing the passive static subsystem, we obtain the static linear 
transformation H(S) (£ ou t = H(S)£) of the dual-NOPA coherent feedback static system as 
a function of S', 


H(S) = H 2 + h 1 P[a?SH 1 H BS ], 


(25) 
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where 


P 


Hbs 

Hi 

h 2 

hi 

h 2 


{h-a 2 S{I 2 ®h 2 ))-\ 

(3 [ 04x2 I A 04x2 ] + Oi(3S 


1 0 0 0 0 0 0 0 

0 1 0 0 0 0 0 0 

0 0 0 0 0 0 1 0 ’ 

0 0 0 0 0 0 0 1 


h\ 0 2 h 3 hi 0 2 0 2 
0 2 hi 0 2 0 2 hi h 3 


h 2 

°2 

hi 

h 3 

o 2 

°2 ' 

04x8 

. O 2 

h 2 

o 2 

o 2 

hs 

hi 


I 2 < 8 > hi, 
h 2 0 
0 -h 2 ’ 



1 2 <E> f ki , 

hi 0 
0 —hi 


(26) 


4 Optimization of S 


In this section, we aim to optimize the EPR entanglement generated in by the dual- 
NOPA coherent feedback system of Fig. [4j in the infinite bandwidth limit, by finding a 
complex unitary matrix at which the two-mode squeezing spectra of the two outgoing 
fields are locally minimized, with respect of S. Since the system is infinite bandwidth, 
V±{mS) = V±(0) for all uj, thus we shall denote V(mj) and V±(iuj) simply as V and V±, 


with no dependence on c 0 . Based on 
spectra is 


and (flTj) , the sum of the two-mode squeezing 


where 


V = V+ + V- 

= Tr [H*Hi + H*H 2 \, 
= Tr [H{S)*M lt2 H{S)] 

0 10 ' 
1 0 -1 
0 1 0 
0-101 



(27) 


(28) 


As V is a function of S or S, we define V(S) as the value of V for a fixed value of S, and 
V(S) as the value of V for a fixed value of S. 

We aim to find a complex unitary matrix S as a local minimizer of the cost function 
V(S). The optimization problem with a unitary constraint can be solved by the method 
of modified steepest descent on a Sticfcl manifold introduced in (23), which employs the 
first-order derivative of the cost function. The Stiefcl manifold in our problem is the set 
St(2, 2) = |S' G C 2x2 : S*S = ^|- 
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Since (/ — Y)~ l = (/ — Y)~ 3 (I + Y — Y) = /+ (/ — Y') -1 F for any square matrix Y such 
that I—Y is invertible, we expand H(S+AS) as H(S)+H(AS)+H(AS 2 )+0(AS 3 ), where 
0(AS 3 ) denotes terms that are products containing at least three AS. H(S), H(AS ) and 
H(AS 2 ) are real matrices, 


H(AS) = PASQ 
H(AS 2 ) = a 2 PAS(L 2 ®h 2 )PASQ, 


where 


Q = a 2 hi 


(it + a‘ (h 0 h-A PS '} II, (I 2 0 ii?)PH BS 


(29) 


(30) 


Following (27) and based on the facts that a matrix and its transpose have the same trace, 
we have 


V(S + AS) = Tt[H(S + AS)*M 1>2 H(S + AS)] 

= V(S) + Tr[H (AS)* Mi j2 H (S) + H(S)*M h2 H(AS) + H(AS 2 )*M ly2 H(S ) 
+H(S)*M h2 H(AS 2 ) + H(AS)*M h2 H(AS)} + 0 (|| AS|| 3 ) 

= V(S) + 2 Tt[MH(AS)} + 2 Tr[MH(AS 2 )] 

+ Tr[H(AS)*M li2 H(AS)] + 0(|| AS|| 3 ), (31) 


where M = H(S)*Mi j2 and 0(||AS|| 3 ) denotes that the function 0(||AS'|| 3 ) satisfies 
°ll^ii* — c f° r some positive constant c for all HASH > 0 sufficiently small. Further¬ 
more, based on (23), we obtain that 


V(S + AS) 


V(S) + ReTi[AS*Dg] + i 


vec(AS) 

* 

x 

vec(AS) 

_ vec(AS # ) _ 

A 

_ vec(AS # ) _ 


+ 0(||AS|| 3 ),(32) 


where 

Dg = 2 K*(QMP) t K, 

x = i [ (R* ® k) ( k * ® k)* ] * h [ (k* ® k) (k# ®k)*], 

h = WL T (QMP) T ®((I 2 ®h 2 )P) + 2(QQ T )® (P t M 1 j 2 P), (33) 
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0 0 
0 0 
1 0 
0 0 
0 0 
0 0 
0 1 
0 0 


0 0 
0 0 
0 0 
1 0 
0 0 
0 0 
0 0 
0 1 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
1 0 
0 0 
0 0 
0 0 
0 1 


Dg is the directional derivative of V(S ) at S in the direction AS [2j]• 


Theorem 4.1 The matrix S c fb corresponding to the dual-NOPA coherent feedback system 
given by (21) is a critical point of the function V(S). 


Proof. According to (24j . we have a one-to-one corresponding cost function g§(AS) on 
the tangent space to the Stiefel manifold at the point S, with AS a vector on this tangent 
space, defined by g§(AS) = V{ir(S + AS)), where ir(-) is the projection operator onto the 
manifold. The descent direction Z d in the tangent space at S is 


Z H = SD%S - D*. 


Based on (33), when S = S c fb , D§ becomes 


D c — d o 

&cfb &cfb 


0 1 
1 0 ’ 


(34) 


(35) 


where da is a real coefficient 

£>cfb 

d§ cf» = (1 + 2a 2 x + 2Kx -x 2 + K 2 x 2 ) 3 ^ ^ 

+ (-l + K 2 )x 2 )( 4x(l + 2 Kx + x 2 + K 2 x 2 ) 

+2a 4 x(— 1 + (—1 + K 2 )x 2 ) + a 2 (—1 — 2Kx 

+ 6/1 x 3 + 2K 3 x 3 + x 4 - 2K 2 x 4 + K 4 x 4 ) )). (36) 


Thus, the descent direction Z d is 

Zd,cfb = ScfbDg cfb S c fb ~ Dg cfb = 02 - 


(37) 
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Thus, the gradient of the function gg(AS) at AS = 0 along the tangent space at S c fb is 
grad (WO)) = S c fbD*~ S c f b — Da = 0 2 (see [23][Eq. (27)]), which establishes that SW 

&cfb C J® 

is a critical point. 


Now we check the Hessian matrix of the function gg(AS). Based on Proposition 12 in 
and (32), we have following the second order expansion along any direction A S on the 
tangent space at S, 


9s(^) 


V(tt(S + AS)) 

V(S) +ReTr[AS*Dg] + ^ 

+0(||A^|| 3 ), 


vec(AS) 

vec(AS' # ) 


Hess (5) 


vec(AS') 

vec(AS ,# ) 


(38) 


where 


Hess(S') = X — - 


(S*Dg) T <g> I 2 0 4 

0 4 ({S*Ds) T ®I 2 )* 


(39) 


denotes the Hessian matrix of gg(AS). Firstly, we consider the system in an ideal case, 
where there are no losses (k — 0 and a — 1). As reported in [6], in this lossless scenario 
the range of x over which the dual-NOPA coherent feedback system is stable in the finite 
bandwidth case is x G [0, a/2 — 1), independently of the actual bandwidth of the NOPAs. 
Thus, it is natural to also take this as the range of admissible values for x in the infinite 
bandwidth limit of this paper. By checking eigenvalues of the Hessian matrix, we have the 
following theorem. 


Theorem 4.2 In the absence of transmission and amplification losses, S c fb is a local min- 
imizer of the function V(S) when x G (\/5 — 2, a/2 — 1). 

Proof. Let a = 1 and k = 0. With the help of Mathematica, the eigenvalues of Hess(S') 
at S — S c fb can be found to be 

8 a; (1 — a/)(l + x 2 ) 2 

61 (1 + 2x — x 2 ) 4 
8x(l — a: 2 )(l + x 2 ) 2 

62 = (1 - 6a; 2 + x 4 ) 2 ’ 

8 a; (1 + x 2 ) 2 (— 1 + 4a; + x 2 ) 

63 (1 + 2a; — x 2 ) 4 

8 a; (1 + a; 2 ) 2 (3 — 6a; + 2a; 2 + 6a; 3 + 3a; 4 ) . . 

64 = (1 + 2a; — a; 2 ) 3 (l + 2a; + a; 2 ) 2 ' ( ) 

As x G (0, \pl — 1), e 1} e 2 and e 4 have positive values, while e 3 > 0 when — 1 + 4x + x 2 > 0, 
that is, x > y/5 — 2. Therefore, for x G (\/5 — 2, \/2 — 1), the Hessian matrix Hess(S^ c /;,) is 
positive definite, which establishes that S c fb is a local minimizer for these values of x. 
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Table [l] and Table [2] illustrate the effect of transmission and amplification losses on 
the range (xi m , a/2 — 1) of over which S c fb is a local minimizer. We see that as either 
transmission losses or amplification losses increase, the range of values of x over which the 
dual-NOPA coherent feedback network is optimal become wider. 


Table 1: Influence of transmission losses on the range (xi m , \/2 — 1) of over which S c fb is 
a local minimizer with k = 0 and a = io~°- 005d _ 


d 

%Inn 

0 

0.236068 

1 

0.212692 

5 

0.134477 


Table 2: Influence of amplification losses on the range (xi m , a/2 — 1) of over which S c fb is 
a local minimizer with d — 1 and a = io~°- 005d _ 


K 

%lm 

0 

0.212692 

n i 3xlO b „ 
U -\/2x0.6 X 

0.211836 

n n 3xlO b 

U -^/2x0.6 X 

0.210989 

V2x0.6 

0.208503 

3x 10 H 

P2x0.6 X 

0.204528 


5 Conclusion 

This paper has studied the optimization of EPR entanglement of a static linear quantum 
system that is composed of a static linear passive optical network in a certain coherent 
feedback configuration with two NOPAs in the infinite bandwidth limit. We reformulate 
the optimization of the EPR entanglement to the problem of finding a 2 x 2 complex 
unitary matrix at which a cost function V(S) is locally minimized, with respect of S. By 
employing the modified steepest descent on Sticfcl manifold method, we have found the 
unitary matrix S c fb corresponding to the coherent feedback system shown as Fig. [i] as a 
critical point of V(S). When losses are neglected, the coherent feedback system is a local 
minimizer when x G (Vb — 2, \/2 —1). When transmission and amplification losses increase, 
the range of values of x over which the coherent feedback system is a local minimizer of 
R(s) is enlarged. In addition, one may wonder if there exists other local minimizers at 
which the system generates better EPR entanglement. Hence future work can consider 
further developing the static passive optical network to search for another local optimizer 
that may yield better EPR entanglement than the system studied in [6] as shown in Fig. [2j 
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